NAR Molecular Medicine

We present a novel integrative analysis of transposable elements (TEs) in 4 single cell RNA-seq (scRNA-seq) datasets of postmortem substantia nigra pars compacta samples of Parkinson Disease (PD) patients matched healthy controls, with the objective of building a cell-type specific trustworthy atlas of TEs that may clarify the role of TEs in sex differences in PD. We have used the soloTE tool to evaluate the TEs expression changes across all snRNA-seq studies identified in our previous systematic review, and then integrated the results using meta-analysis techniques. Finally, we evaluated the possible associations between TEs and protein coding genes by integrating our previous results in this matter with the information of TEs obtained, in order to propose the possible action mechanism by which some of the TEs contribute to PD.

The role of the epigenome in age-related neurodegenerative disorders remains understudied. Here, we analyzed circulating cell-free DNA (cfDNA) from blood to detect methylation changes as a liquid-biopsy for Amyotrophic Lateral Sclerosis (ALS). Our study included 20 patients with sporadic ALS, 10 patients with C9orf72-associated ALS, 10 asymptomatic carriers of the C9orf72 repeat expansion mutation, and 21 non-disease controls. Following targeted enzymatic methyl-sequencing (EM-seq) of [~]4 million CpG sites, we detected numerous differentially methylated genes, including several implicated in ALS disease risk and pathogenesis. By integrating multiple epigenetic features, we delineated a distinct epigenetic signature, which achieved an average area under the curve (AUC) of 0.91 {+/-} 0.10 upon receiver operator characteristic (ROC) analysis, which enabled detection of [~]70% of ALS patients with close to 100% specificity. Furthermore, we also identified a set of genes whose methylation status significantly correlated with clinical disease progression and cerebrospinal fluid (CSF) neurofilament levels. Our results reveal the potential of cfDNA-based biomarkers to accurately diagnose ALS and potentially predict disease progression.

Long-read RNA sequencing (lrRNA-seq) provides advantages for transcript discovery and quantification through the sequencing of full-length transcripts. Although recent benchmarks have evaluated long-read technologies and quantification tools, to the best of our knowledge, no study to date has jointly compared sequencing technology, quantification choice, and depth across both bulk and single-cell platforms. Here, we generate a matched dataset using NGN2-induced neurons derived from Fragile X syndrome and isogenic rescue lines, profiled with bulk and single-cell Illumina, Oxford Nanopore Technologies (ONT), and Pacific Biosciences (PB) Kinnex technologies. All platforms and technologies capture the expected FMR1 reactivation signal. We find that PB bulk under-detects and under-quantifies short transcripts (less than 1.25 kb), ONT bulk under-detects and under-quantifies long transcripts (greater than 5 kb), and single-cell long-read technologies a large number of single-cell specific transcripts associated with truncations. Across six bulk and four single-cell long-read quantification tools, Isosceles, Miniquant, and Oarfish provide the best compromise between computational efficiency and performance in terms of quantification accuracy as measured by spike-ins, comparisons to Illumina, and on consensus based down-stream tasks such as differential transcript expression (DTE). Depth-equivalency analyses reveal that PB single-cell sequencing requires approximately three- to four-fold greater depth than bulk to reach comparable power for transcript discovery and differential transcript expression. Our work aims to offer practical guidance for study design, including the choice of technology, sequencing depth, and quantification method. In addition, we hope our data may serve a reference dataset to evaluate emerging long-read transcriptomic protocols and methods as well as more closely investigate FMR1 biology.

Accurate detection of RNA splice variants is often hindered when transcripts lack large distinguishable exonic regions, making conventional PCR strategies challenging. We developed a simple melting temperature (Tm)-guided exon-exon junction (EEJ) RT-PCR method to enable variant-specific detection under these conditions. Uni-directional primers spanning exon-exon junctions were designed so that approximately each half anneals to adjacent exons. The Tm of each half-site was set >7{degrees}C below the annealing temperature, preventing stable binding to individual exons and enforcing junction-dependent amplification. The method was evaluated using HTRA1-AS1 long noncoding RNA variants that share overlapping exon sequences but differ in splice connectivity. HTRA1-AS1 comprises five variants, only one with a large distinguishable exon. Tm-guided EEJ primers robustly discriminated the remaining four variants. After optimization, amplification yielded sharp, single bands with minimal cross-reactivity. Compared with conventional designs, this approach reduced heteroduplex and heteroquadruplex formation, improving band clarity. Sanger sequencing confirmed junction specificity, and the method performed well in multiplex settings. Overall, Tm-guided EEJ RT-PCR is a cost-effective, high-resolution approach for detecting RNA variants lacking easily distinguishable exonic regions, readily compatible with standard RT-PCR and qPCR workflows.

The CRISPR-Cas9 system has revolutionized genome engineering, yet its full therapeutic potential remains constrained by challenges in precisely modulating its activity and specificity. Here we report a fully computational end-to-end pipeline for the de novo design of a single-domain VHH nanobody (NbSpCas9-v1) targeting a structurally conserved, non-catalytic epitope at the PAM-interacting (PI) and RuvC-III interface of Streptococcus pyogenes Cas9 (SpCas9; PDB: 4UN3). Nanobody sequences were generated using BoltzGen, a generative diffusion binder design framework, and co-folded with SpCas9 using Boltz-2 to evaluate structural confidence and binding affinity. The top-ranked model (SpCas9_4UN3_Bivalent_Hub_v1) achieved a complex pLDDT of 0.8406, an aggregate score of 0.8016, and an ipTM of >0.8, indicating high confidence in the nanobody-antigen interface. The designed 1,616-residue quaternary complex (SpCas9 + sgRNA + DNA + nanobody) was subjected to 10 ns of all-atom molecular dynamics (MD) simulation using the AMBER14SB force field within the GROMACS/OpenMM framework. The complex stabilized at RMSD [~]6 [A] with a radius of gyration of 39-44 [A], confirming thermodynamic stability under physiological conditions (310 K, 0.15 M NaCl). A conserved 96.3 [A] inter-molecular distance between the nanobody centroid and the HNH catalytic residue H840 establishes NbSpCas9-v1 as a distal, non-inhibitory binder -- ideally suited for a Bivalent Hub architecture recruiting secondary effectors to the Cas9 ribonucleoprotein (RNP). The nanobody-Cas9 interface is stabilized by 8 hydrogen bonds, 4 salt bridges, and [~]1,850 [A]2 of buried solvent-accessible surface area. These results provide a rigorous structural and dynamic foundation for experimental validation of VHH-based CRISPR-Cas9 enhancers and modulators. GRAPHICAL ABSTRACTThe computational workflow proceeds from SpCas9 crystal structure acquisition (PDB: 4UN3) through BoltzGen nanobody design, Boltz-2 structural co-folding, 10 ns explicit-solvent MD validation, and Bivalent Hub functional characterization. The PyMOL rendering below shows the full quaternary complex at atomistic resolution.

Global quantification of DNA cytosine modifications, including 5-methylcytosine (5-mC) and 5-hydroxymethylcytosine (5-hmC), is important for understanding cancer biology, though established methods require multi-step workflows and costly instrumentation. Here we show that attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy combined with regression modelling enables rapid, label-free, and non-destructive quantification of both modifications from DNA samples. Using Adenomatous Polyposis Coli (APC) promoter DNA standards spanning 0-100% modification, we identified modification-sensitive spectral features and observed that 5-hmC produces greater spectral changes than 5-mC. A univariate peak-ratio approach yielded strong linearity for both modifications (R2 = 0.97), while partial least squares regression (PLSR) improved quantification accuracy to R2 = 0.99 (RMSE = 2.6%) for 5-hmC and R2 = 0.97 (RMSE = 5.7%) for 5-mC. In composite mixtures containing all three cytosine states, 5-hmC remained highly quantifiable (R2 = 0.97; RMSE = 5.1%), while 5-mC accuracy decreased (R2 = 0.90; RMSE = 9.6%), consistent with the greater spectral distinctiveness secondary to the hydroxymethyl group. Transferability was assessed using circulating tumour DNA (ctDNA), short cell-free DNA fragments shed from tumour cells into the bloodstream, comprising multiplexed reference material spanning seven genomic regions and a polydisperse fragment-length distribution (155-220 bp). After domain adaptation between synthetic and ctDNA spectra, we obtained a quantitative methylation calibration with R2 = 0.98 and RMSE = 5.2% under cross-validation. These results support ATR-FTIR spectroscopy as a viable platform for global cytosine modification quantification and establish proof-of-concept applicability to ctDNA analysis.

Spinal muscular atrophy (SMA), traditionally defined as a neuromuscular disorder characterized by degeneration of lower motor neurons, is increasingly recognized as a multi-organ disease. SMA is caused by deficiency of the survival motor neuron (SMN) protein below a critical threshold required for cellular homeostasis. While motor neurons are particularly vulnerable, the ubiquitous expression and fundamental functions of SMN result in widespread perturbations across multiple tissues. Here, we generated a label-free quantitative proteomics atlas of spinal cord, heart, and gastrocnemius muscle from wild-type, heterozygous, and SMA mice at the symptomatic stage, including cohorts treated, at postnatal day 1 (P1), with a systemic suboptimal dose of SMN antisense oligonucleotides (SMN-ASOs), resulting in partial SMN restoration. SMN deficiency induced pronounced, tissue-specific proteome remodeling, with peripheral tissues exhibiting broader molecular alterations than spinal cord. Cross-tissue analyses revealed limited overlap, although heart and muscle showed partial convergence in metabolic and mitochondrial-associated pathways. SMN-ASO treatment partially repositioned these proteomes toward control states; however, restoration was incomplete and strongly tissue-dependent, with persistent dysregulation of mitochondrial and metabolic pathways. These findings demonstrate that SMN deficiency drives systemic yet heterogeneous proteome remodeling and that partial SMN restoration does not fully reverse established molecular alterations. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=99 SRC="FIGDIR/small/715402v1_ufig1.gif" ALT="Figure 1"> View larger version (35K): org.highwire.dtl.DTLVardef@114d73borg.highwire.dtl.DTLVardef@13e8c13org.highwire.dtl.DTLVardef@15e4ba0org.highwire.dtl.DTLVardef@1b70fb8_HPS_FORMAT_FIGEXP M_FIG C_FIG

Mammalian telomeric DNA comprises long tracts of tandem TTAGGG repeats. The same repeats are also found at internal chromosomal regions called interstitial telomeric sequences (ITSs). Telomeres are transcribed into UUAGGG-containing transcripts, named TERRA, which serve multiple functions in maintaining telomere integrity. Complementary RNAs containing C-rich telomeric repeats, named ARIA, have also been identified in few yeast mutants and mammalian cells with dysfunctional telomeres. The molecular features and functions of ARIA remain understudied, mainly due to its low abundance and the lack of suitable cellular systems. Here, we show that Chinese hamster ovary (CHO) cells produce abundant TERRA and ARIA transcripts, predominantly originating from ITSs. Both RNAs are polyadenylated, exhibit relatively short half-lives and form large cellular foci. We also show that ARIA depletion leads to exposure of single-stranded (ss) DNA at ITSs and that ssDNA exposure increases when ITS DNA is damaged. SsDNA formation does not require the DNA damage signaling kinases ATM and ATR, nor the exonucleases DNA2 and EXO1; however, ATM prevents excessive ssDNA accumulation when ARIA function is inhibited. These findings establish CHO cells as a powerful model to dissect telomeric RNA functions and reveal ARIA as a key regulator of telomeric repeat DNA integrity.

Recombinant adeno-associated virus (rAAV) vectors are a leading platform for gene delivery in basic and clinical research, yet large-scale manufacturing remains constrained by residual nucleic-acid impurities that compromise safety. In this study, we profiled the DNA species packaged within rAAV capsids and identified plasmid backbone sequences and host cell genomic DNA (hcDNA) as predominant contaminants. To mitigate this critical quality attribute, we implemented upstream strategies designed to fragment or excise backbone DNA, including TelN/TelROL excision, I-SceI meganuclease digestion, CRISPR/Cas9 cleavage, and Cre/LoxP recombination. Quantitatively, TelN/TelROL and I-SceI reduced encapsidated plasmid backbone DNA to approximately 20-30% and 20-40% of baseline levels, respectively, while CRISPR/Cas9 lowered it to about 10-20%. Notably, the Cre/LoxP system eliminated detectable plasmid backbone DNA without compromising vector-genome titers, indicating preserved genomic integrity. Additionlly, supplementating cell culture with a caspase inhibitor significantly reduced hcDNA contamination in rAAV particles to 1-5% of the baseline level. Collectively, these interventions provide practical bioprocess frameworks that markedly enhance rAAV purity via targeted DNA minimization and prevention of hcDNA fragmentation, thereby strengthening the safety profile of rAAV therapeutics in alignment with current Good Manufacturing Practice (cGMP) expectations.

Human mitochondrial genome (mtDNA) encodes multiple proteins in the oxidative phosphorylation complexes as well as the ribosomal and transfer RNAs (tRNAs) needed for in situ translation. These genes are transcribed from only three promoters, producing polycistronic transcripts that are co-transcriptionally cleaved by mitochondrial RNase enzymes to release majority of individual gene products. tRNAs separate many of these genes and are thought to serve as "punctuation" marks that enable RNase recognition, binding, and hydrolysis of the 5' "leader" and 3' "trailer" sequences flanking the tRNA. Mutations in the tRNA genes dominate the mtDNA-linked mitochondrial pathologies; yet a systematic study of the impact of tRNA sequence variation on the RNase-catalyzed processing is lacking. Here, we employed human mitochondrial tRNATyr as a model system to dissect the effect of tRNA variants on the in vitro 5' leader and 3' trailer hydrolysis. We found that nucleotide variations located near the catalytic interfaces - particularly within or near the tRNA acceptor stem - showed the strongest defects in 5' processing and prevented release of the downstream tRNA in a tRNA cluster where multiple tRNAs are transcribed in tandem. This work provides mechanistic insight into how mutations disrupt coordinated mitochondrial tRNA processing and establish a framework for predicting variant effects based on their structural position relative to the processing enzymes.

Alzheimers disease (AD), the leading cause of dementia, affects over 33 million people worldwide, with pathogenesis tied to amyloid-{beta} (A{beta}) accumulation. Although anti-A{beta} monoclonal antibodies have shown clinical benefits, they often cause side effects including amyloid-related imaging abnormalities and brain microhemorrhage, especially in APOE E4 allele carriers. Here we used PHP.eB serotype adeno-associated virus (AAV), a vector with enhanced central nervous system (CNS) tropism, to deliver an A{beta} antibody expression vector (AAV-LEC) into the CNS of APP/PS1 and 5xFAD mice intravenously. The AAV-LEC-mediated expression of anti-A{beta} antibodies in the CNS significantly reduced the number and size of A{beta} plaques at various stages in both APP/PS1 and 5xFAD mice, alongside improved spatial learning and memory. It also reversed abnormal glial activation with reduced disease-associated microglia and astrocytes, and restored oligodendrocyte differentiation and myelin formation. No brain microhemorrhage or liver damage was detected following the AAV-antibody treatment. Thus, this AAV-mediated strategy offers a promising, convenient and safe AD therapeutic approach in the future.

Gene regulation through translation is critical for spatiotemporal protein expression. Internal ribosomal entry sites (IRESes) mediate mRNA-specific translation by recruiting ribosomes to 5 untranslated regions. Circular RNAs (circRNAs), naturally occurring and stable RNA species, are increasingly used as synthetic tools for sustained therapeutic protein translation by IRES-driven initiation. However, the functionality of different IRESes in synthetic circRNAs remains sparsely characterized. We systematically examine circRNA reporter translation by viral and cellular IRESes in human cells and in diverse in vitro translation systems. Improved circRNA purification by urea-PAGE and RNase R-treatment removes contaminants that induce RNA sensing. Viral CVB3 and HCV, as well as cellular Hoxa9, Chrdl1, Cofilin and c-Myc IRESes, effectively drive circRNA translation. We also establish circRNA translation in an improved human cell-free extract that recapitulates IRES-dependent regulation, and allows for precise engineering of HCV IRES-mediated translation. These findings inform IRES selection for synthetic circRNA translation relevant for circRNA-based medicines.

Alzheimers disease (AD) is the most prevalent form of dementia, characterized by progressive memory loss, cognitive decline, and emotional dysregulation. Adult hippocampal neurogenesis (AHN) critically contributes to cognition and mood but undergoes precipitous decline during AD progression. Here, we investigated whether enhancing AHN through genetic expansion of endogenous neural stem cells (NSC) ameliorates AD-related phenotypes. Using lentiviral overexpression of the cell cycle regulators Cdk4 and CyclinD1 in the dentate gyrus of 3xTg-AD mouse, we show that enhancing AHN partially rescues hippocampal-specific cognitive functions, namely: spatial navigation and exploratory behavior. These findings show that endogenous NSC can be exploited to ameliorate hippocampal cognitive functions in AD, providing additional evidence for exploiting AHN as a promising therapeutic target for neurodegenerative disease.

Rotavirus is a major cause of severe diarrheal disease in children under the age of five, with reduced vaccine effectiveness in low-resource settings causing substantial morbidity and mortality. In the absence of approved antiviral therapeutics, treatment is largely supportive, urging the need for targeted and precision-based interventions. VP4 protein plays an essential role in viral attachment, entry, and infectivity, making it a suitable target for targeted therapy. In this context, RNA interference is a specific method for inhibiting viral gene expression with its efficacy depending on sequence conservation, target accessibility, and compatibility with the RISC-loading machinery. In the present study, an integrative in silico approach was employed to design and evaluate siRNAs targeting conserved regions of the VP4 gene across six geographically diverse countries. Candidate siRNAs were screened using established design rules and regression-based scoring with off-target filtering. Three optimized siRNAs were further assessed through structural modeling, molecular docking, and molecular dynamics simulations to examine interactions with human Dicer, TRBP, and Argonaute-2. Comparative dynamic analyses identified one siRNA with enhanced structural compatibility, reduced conformational fluctuations, and stable interactions with RISC-loading proteins. These findings provide a rational computational basis for VP4-targeted siRNA development, facilitating experimental validation.

Muscleblind-like (MBNL) RNA-binding proteins (RBPs) possess modular domains that mediate regulation of alternative splicing and RNA localization. Myotonic Dystrophy Type 1 is a CTG repeat expansion disorder where MBNL is sequestered into intranuclear RNA foci, impairing its function. Previous studies found that MBNL self-associates through its exon 7, but the nature of this interaction is not well understood. We identified a cysteine in MBNL1 exon 7 that enables dimerization through formation of an intermolecular disulfide bond. We likewise demonstrate that MBNL2 dimerizes by forming disulfide bonds between multiple cysteines in its carboxy-terminus. Nucleocytoplasmic fractionation revealed a greater proportion of MBNL1 dimer in the nucleus, suggesting a nuclear function for the MBNL1 dimer. We investigated a connection between MBNL1 dimerization and MBNL1-mediated regulation of alternative splicing. To accomplish this, we mutated the MBNL1 cysteine in question to alanine (C325A) and performed RNAseq. We uncovered novel splicing events sensitive to MBNL1 dimerization. We also found that MBNL1 C325A, when co-expressed with expanded CTG repeats, produces smaller, more numerous foci, suggesting a role for the MBNL1 dimer in maintaining foci integrity. These results provide insight into biological and pathological mechanisms of MBNL1 dimerization and suggest other RBPs might similarly dimerize to regulate function. GRAPHICAL ABSTRACT

Background: Mucopolysaccharidosis type IIIA (MPS IIIA; Sanfilippo syndrome) is a fatal neurodegenerative lysosomal storage disorder caused by impaired degradation of heparan sulfate (HS). Despite rapid advances in gene and enzyme therapies, there remains a critical need for an analytically validated, quantitative biomarker that accurately reflects central nervous system (CNS) substrate burden. Such biomarker would be a valuable tool in assessing disease progression and monitoring therapeutic efficacy. Objective: This study describes the method development, fit for purpose validation, and preliminary clinical application of a quantitative liquid chromatography-mass spectrometry (LC-MS/MS) assay for the HS-derived disaccharide N-sulfoglucosamine-glucuronic acid (GlcNS-GlcUA) in human cerebrospinal fluid (CSF), a critical biomarker for diagnosis, disease monitoring, and regulatory evaluation of emerging MPS IIIA therapies. Methods: A structurally defined GlcNS-GlcUA reference standard and its [13C6]-labeled internal standard were used in a derivatization and detection workflow employing 1-phenyl-3-methyl-5-pyrazolone labeling, and LC-MS/MS. Results: The method exhibited acceptable linearity across 0.005-0.500 nmol/mL (r[≥]0.9976), with intra- and inter-assay imprecision [≤]3.5%CV and accuracy within 95%-110% of nominal concentrations. No matrix or hemolysis interference or carryover was observed, and the analyte remained stable during freeze-thaw storage conditions. Application of the method to 12 CSF samples from patients with MPS IIIA demonstrated quantifiable GlcNS-GlcUA levels ranging from 0.0054 to 0.106 nmol/mL, confirming suitability for clinical and regulatory use. Comparison of the MPS IIIA sample results between the development laboratory and the contract research organization laboratory support robust inter-lab assay transfer. Conclusions: This validated LC-MS/MS method establishes a regulatory-grade quantitative assay for measurement of CSF HS in MPS IIIA. Its high analytical sensitivity and reproducibility enable reliable assessment of CNS substrate reduction and pharmacodynamic response, supporting biomarker-driven therapeutic development and accelerated approval pathways for neuronopathic mucopolysaccharidoses.

Single molecule methods have become prevalent tools in elucidating molecular processes across various life science fields over the past three decades, driving breakthroughs in understanding their underlying molecular mechanisms. In our study, we employed two single-molecule methods, Forster Resonance Energy Transfer (smFRET) and high-resolution optical tweezers, to investigate the transcription of Mycobacterium tuberculosis RNA polymerase (MtbRNAP) from initiation through to termination. We aim to provide a set of comprehensive biophysical tools to deepen our current understanding of MtbRNAP and its transcription factors. These experimental assays represent an important step towards unraveling the molecular dynamics and interactions that support transcription in Mycobacterium tuberculosis.

Spatial transcriptomics (ST) has emerged as a transformative tool for resolving the molecular heterogeneity of complex tissues within their native anatomical context. However, next-generation sequencing (NGS)-based ST platforms frequently encounter sensitivity bottlenecks arising from sub-optimal probe architectures on solid substrates. Conventional single-stranded DNA coupling methods often lead to disordered interfacial molecular conformations due to non-specific nucleobase-mediated surface tethering, which creates steric hindrance and inhibits the enzymatic efficiency of in situ library preparation. Here, we present Decomap (double-strand protected combinatorial barcoding microarray chip), a high-performance ST platform utilizing a triple-segment (dsZ-X-Y) fabrication strategy to achieve superior transcript capture efficiency. This structural optimization significantly enhances DNA ligation kinetics and subsequent polymerase-mediated extension, overcoming the fundamental limitations of traditional single-stranded coupling strategies. Decomap-seq achieved a median detection of 7,200 genes and 29,097 UMIs per 50 m-spot at a sequencing saturation of 50.1%. These results validate Decomap as a highly sensitive and robust tool for spatial transcriptomics, offering a powerful platform for advancing research in histopathology, developmental biology, and neuroscience.

TDP-43 is a nuclear RNA-binding protein regulating numerous steps in RNA metabolism, including alternative splicing. It is a major pathological hallmark of several neurodegenerative diseases, where it forms cytoplasmic aggregates in affected brain regions. TDP-43 can undergo phase separation (PS) and this condensation behavior may be linked to aggregate formation. Whether and how PS governs TDP-43s RNA regulatory functions remains poorly understood. Here we utilized rationally designed mutations in the TDP-43 low complexity domain to tune TDP-43 PS, yielding a panel of TDP-43 variants with reduced propensity to form condensates ("PS-deficient"), and a panel of TDP-43 variants that form more irreversible, undynamic condensates ("solid-like") in vitro and in cells. Affinity proteomics coupled to mass spectrometry identified a set of interactors whose association with TDP-43 is PS-dependent. This includes multiple splicing regulators and the RNA helicase UPF1, which show increased interactions with solid-like TDP-43 variants. Consistently, we identified PS-dependent alternative splicing events that translate into measurable changes in RNA and protein abundance. Our results highlight that TDP-43 PS regulates RNA and protein homeostasis both directly, by altering a subset of TDP-43-dependent alternative splicing events, and indirectly, by changing interactions with other RNA regulatory factors.

Transposable elements (TEs) are mobile genetic sequences that can generate new copies of themselves via insertional mutations. These viral-like sequences comprise nearly half the human genome and are present in most genome wide sequencing assays. While only a small fraction of genomic TEs have retained their ability to transpose, TE sequences are often transcribed from their own promoters or as part of larger gene transcripts. Accurately assessing TE expression from each individual genomic TE locus remains an open problem in the field, due to the highly repetitive nature of these multi-copy sequences. These issues are compounded in single-cell and single-nucleus transcriptome experiments, where additional complications arise due to sparse read coverage and unprocessed mRNA introns. Here we present our tool for single-cell TE and gene expression analysis, TEsingle. Using synthetic datasets, we show the problems that arise when not properly accounting for intron retention events, failing to address uncertainty in alignment scoring, and failing to make use of unique molecular identifiers for transcript resolution. Addressing these challenges has enabled an accurate TE analysis suite that simultaneously tracks gene expression as well as locus-specific resolution of expressed TEs. We showcase the performance of TEsingle using single-nucleus profiles from substantia nigra (SN) tissues of Parkinsons Disease (PD) patients. We find examples of young and intact TEs that mark dopaminergic neurons (DA) as well as many young TEs from the LINE and ERV families that are elevated in PD neurons and glia. These results demonstrate that TE expression is highly cell-type and cellular-state specific and elevated in particular subsets of neurons, astrocytes, and microglia from PD patients.